Quantum Spin Hall Effect and Topologically Invariant Chern Numbers 
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We present a topological description of quantum spin Hall effect (QSHE) in a two-dimensional 
electron system on honeycomb lattice with both intrinsic and Rashba spin-orbit couplings. We show 
that the topology of the band insulator can be characterized by a 2 x 2 traceless matrix of first Chern 
integers. The nontrivial QSHE phase is identified by the nonzero diagonal matrix elements of the 
Chern number matrix (CNM). A spin Chern number is derived from the CNM, which is conserved 
in the presence of finite disorder scattering and spin nonconserving Rashba coupling. By using the 
Laughlin's gedanken experiment, we numerically calculate the spin polarization and spin transfer 
rate of the conducting edge states, and determine a phase diagram for the QSHE. 

PACS numbers: 73.43.-f, 72.25.-b, 73.43.Nq 



Topological quantities are fundamentally important in 
characterizing a new kind of order — topological ordeni, 
which cannot be described in the paradigm of a Landau 
symmetry breaking theory. For two-dimensional (2D) 
electron systems in perpendicular magnetic fields, they 
are manifested in the transverse electrical transport in in- 
teger and fractional quantum Hall effect (IQHE/FQHE) 
states^. It was first revealed by Thouless et. &l£ that 
each IQHE state is associated with a topologically invari- 
ant integer known as first Chern number, which precisely 
equals the Hall conductance in units of e 2 /h. The exact 
quantization of the Hall conductance can also be formu- 
lated in terms of a 2D band-structure Berry phase^ii^, 
which remains an integral invariant till the band energy 
gap (or more precisely the mobility gap£ in the presence 
of disorder) collapses. 

While the conventional IQHE is usually associated 
with strong magnetic fields, Haldane^ has explicitly 
shown that it can actually occur in the absence of 
magnetic field in band insulators with graphcnc-likc 
band structure. The one-component Haldane's model 
explicitly breaks time-reversal symmetry, resulting in 
a condensed-matter realization of a parity symmetry 
anomaly with chiral edge states at the boundary of the 
sample. In realistic electron systems, however, the cou- 
pled spin degrees of freedom can recover the time reversal 
symmetry by forming Kramers degenerate states, which 
belong to the universality class of zero charge Chern num- 
ber as the total Berry curvature of the occupied energy 
band of both spins sums to zero. 

This class of insulator has been recently found 9 ^ 
to possess a dissipationlcss quantum spin Hall effect 
(QSHE)ii, which is distinct from the intrinsic spin Hall 
effect (SHE) in the metallic systems^. The QSHE has 
been shown to be robust against disorder scattering and 
other perturbation effect a 9 ' 10 . Whether there exists an 
underlying topological invariant "protecting" the QSHE 
is a very important issue for both fundamental under- 



standing and potential applications of the QSHE. While 
the previously proposed 9 *^ classification of the QSHE 
suggests that the conducting edge states are protected by 
time-reversal symmetry, it does not distinguish between 
two QSHE states with spin Hall conductance (SHC) of 
opposite signs. Thus it remains an open issue if the 
QSHE states can be classified by a more definite topo- 
logical quantum number similar to the Chern number for 
the conventional IQHE. 

In this Letter, we clarify the topological nature of the 
QSHE in graphene model with spin-orbit coupling. We 
present numerical evidence that a pair of nonchiral con- 
ducting edge states are responsible for the dissipationlcss 
spin current in the open system with disorder scatter- 
ing. The spin transfer and spin accumulation are calcu- 
lated, which increase linearly with the flux insertion in 
the Laughlin's gedanken experiment^. We further es- 
tablish the relation between the conducting edge states, 
the topological invariant Chern number matrix (CNM) 
and a new conserved spin Chern number for the bulk 
system with electron Fermi energy lying inside the band 
gap. The nontrivial spin Chern number distinguishes a 
QSHE state from an ordinary band insulating state, and 
is responsible for the robust QSHC. 

We begin with the 2D honeycomb lattice mode&2*iS, 
which is relevant to the 2D electrons in a single-atomic 
layer graphene system^: 

+ iVR 4% ' ( a X d 'j) C 3 + X! WiC i C i i I 1 ) 
(ij) » 

where cj = (cL,c|.) are electron creation operators, and 
er is the Pauli matrix. The first term is the usual near- 
est neighbor hopping term, and the second term is an 
intrinsic SO coupling preserving the lattice symmetries 
with i and j as two next nearest neighbor sites and k 
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FIG. 1: The evolution of edge levels inside the bulk gap with 
twist boundary phase 9 y in a cylindric 2D sample of N = 
480 x 240, Vso = 0.05t, Vr = 0.05t and W = l.Ot. 

their unique common nearest neighbor. Here, the vector 
djfc points from k to i, with the distance between two 
nearest neighbor sites taken to be unity. The third term 
stands for the Rashba SO coupling with strength Vr, and 
the last term represents an on-site random disorder with 
strength \w z \ < W/2. 

In the absence of the Rashba SO coupling Vr = 0, 
the model Eq. Q reduces to a two-component Haldane's 
modei&, which exhibits a quantized SHC <j s h = ±2^- 
with sign depending on the sign of Vso- Such a QSHE 
has been shown to persist into a finite range of Vr~ and 
disorder strength W— , so long as the bulk band gap re- 
mains finite. In the following, we will perform systematic 
numerical simulations designed to reveal the topological 
characteristics of such a nontrivial band insulator. 

Gauge experiment and edge states free from backward 
scattering — Let us first study a rectangular geometry of 
honeycomb lattice consisting of N y zigzag chains with N x 
atom sites on each chain. A twisted boundary condition 
is imposed^ in the y-direction, and an open boundary 
condition is used in the ^-direction. The system is also 
equivalent to a cylindric geometry threaded by a flux <3? = 
in the axial direction (i-dircction), similar to that 
used in the Laughlin's gedanken experiment— for IQHE. 
The calculated energy levels for Vso = 0.054, Vr = 0.054, 
W = l.Ot and system size N = N x x N y = 480 x 240 are 
plotted in Fig. ^ We see that some new states emerge 
within the bulk energy gap — O.lf < E m < 0.2t with 
the density of states much smaller than that of the bulk 
states outside the energy gap. These mid-gap states can 
be identified as edge states with their wavefunctions lo- 
calized within a few lattice constants at the two open 
boundaries in the i-direction and extended along the y- 
direction. 

The edge states in Fig. ^ continually evolve and cross 
each other with increasing 8 y . Similar to IQHE, the level 
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FIG. 2: Spin polarization P„ defined in the text (in units of 
h/2) as a function of eigen energy E m for two random disorder 
configurations (we have checked over 100 configurations) of 
strength W = l.Ot at different system sizes N = 120 x 120 — 



480 x 480 (with N x 
in Fig. 1. 



N y = L). Vso and Vr are the same as 



crossing of the edge states in Fig. ^ clearly demonstrates 
the absence of backward scattering, and thus each state 
will not simply evolve back to its original state after the 
insertion of one flux quantum $ = 1 (or 9 y = 2ir) , which 
leads to an adiabatic transverse transport of electrons 
from one edge to the other—. For a Fermi energy Ef 
inside the band gap, as indicated by the dashed line in 
Fig. there are always a pair of edge states crossing the 
Ef line and evolving out of the Fermi sea. They can be 
identified in the real space as moving oppositely, cither 
from the left edge to the right one or vice versa. Thus 
these ideal conducting channels arc nonchiral, unlike the 
IQHE, and do not contribute to a net transverse electric 
charge transport. 

Robust spin polarization carried by the edge states — 
We then explicitly compute the spin polarization = 
f ( TO I Si cfa z Ci\m) along the z-direction with \m) as 
the m-th single-particle eigenstate. Spin polarizations 
along the x- and y-directions are found much smaller. 
At y = 0, all the eigenstates form Kramers degener- 
ate pairs, and each pair carry totally zero spin due to 
time-reversal symmetry. We find that upon insertion of 
the flux (6 y 7^ 0), a nonzero spin polarization is quickly 
developed for each edge state. Fig. |21 shows the spin po- 
larization of each eigenstate as a function of cigenencrgy 
E m at 9 y = 27r/48. The spin polarization is finite only 
for the edge states within the energy gap and insensitive 
to disorder configurations and sample sizes, as all the 
data points nicely collapse together. We find that is 
mainly affected by the Rashba coupling Vr and decreases 
monotonically with increasing Vr. 

Spin transfer rate and phase diagram — Now we can 
determine the spin transfer by the edge states with the 
adiabatic insertion of the magnetic flux. There are always 
a pair of states crossing the Ef = line from below, as 
shown in Fig.Q] One state carries a positive P^, which 
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FIG. 3: (a) Spin transfer rate 7 = A ^ 5 A ^ dgc (in units of 
e/4-7r) vs. disorder strength W at Vso = 0.05f and Vr = 
0.05t averaged over 500 disorder configurations; (b) The phase 
boundary for the QSHE determined from the spin transfer, 
where the transition happens with a change of the numbers 
of edge channels from N e d = 2 to N e d = 0. 

moves to the right edge of the sample, and the other with 
negative moves to the left edge, as determined from 
the wavefunctions. Thus after the insertion of one flux 



quantum 



— ► 2tt), a net spin transfer occurs, which 



results in spin accumulations at the two open edges of the 
sample. We can determine the average spin transfer rate 

7 by using the relation 7 = A £ ds ° = J d6 y d g° ds ° 

with (S' 2 ) dgo = Yl'm^rm where the summation runs 
over only the edge states at the right edge. As shown 
in Fig. [3^, the spin transfer rate remains about 1.8 (in 
units of j-) till a critical disorder strength W c ~ 2 At 
is reached, beyond which 7 vanishes quickly. By anal- 
ogy with IQHE—, the spin transfer rate corresponds to a 
finite SHC of value around 1-8^, which is in good agree- 
ment with the value obtained by using the Landauer- 
Biittiker formulaiS. 

The robust and dissipationless edge states thus clearly 
distinguishes the QSHE in the present system from the 
intrinsic SHE in other metallic systems. In Fig. [^Jd, a 
phase diagram in the parameter plane of the Rashba cou- 
pling Vr and disorder strength W is obtained, where the 
phase boundary of the QSHE regime is determined by 
numerical calculation of the above spin transfer rate. We 
find that such a phase boundary is well correlated with 
the collapse of the bulk energy gap, although the precise 
position of the latter is harder to pinpoint at finite W. 

Topologically invariant Chern number matrix — At this 
stage, we have firmly established the connection between 
the QSHE in the model Eq. Q and the presence of edge 
states within the bulk band gap. In IQHE, it is well 
known that edge states are closely related to the bulk 
topological quantity, i.e., the first Chern number. In the 
following, we examine whether there also exists a topo- 
logical characterization of the bulk states, which is re- 
sponsible for the existence of the edge states. 

Generally, in diagonalizing the Hamiltonian Eq. 



one may introduce a generalized boundary conditio] 
for a 2D many-body wavefunction: $(..., r.; a + Lj, ...) 
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e l0 i <&(..., Ti a , ...) , where j = x,y, and the system length 
vector Li x = -^-ai and L y = N y a.2, with ai and &2 as 
two primitive vectors of the Bravais lattice. The twisted 
boundary condition is represented by < 6? < 2n, where 
ct =| and I denotes the spin index. Through a unitary 

transformation * = exp[— i Y, a "5Zi a \i^ x i a + z^2/i a )]$! 
where the summation runs over all electrons of both 
spins, 'F becomes periodic on a torus. One can then 
define the topological Chern numbers as&i£: 
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where the area integration is over a unit cell < 9% , 
0§ < 2tt. With a, (3 =t, |, C a ^ form a 2 x 2 CNMii. 

The many-body wavefunction ^(6) is necessarily a 
smooth function of the boundary phase 9 = {9^,9^} 
when the energy gap remains. One can then prove the 
exact quantization of C a ^ by strictly following the ar- 
gument of Thouless et al&. Each topologically invariant 
matrix element C"' 13 should remain unchanged until the 
energy gap collapses. Without breaking the time reversal 
symmetry, the trace of the CNM should always vanish. 
Consider the simple case without Rashba coupling term 
(Vr = 0), the only nonzero matrix elements are the di- 
agonal ones: C a a = ij a with 77-f = — r]i = 1 for Vso > 0, 
which change sign with Vso- They describe the IQHE 
in decoupled subsystems of a =| and J,. While the total 
charge Chern number C c = a C a, @ (corresponding to 
the total charge Hall conductance of the system)i£ can- 
cels out, the total spin-related Chern number is quan- 
tized to C sc = J2 a p a C a 'P = 2, which represents the 
spin Hall response when a common electric field for both 
spin components is imposed along the y-direction. In this 
spin decoupled limit, the quantized spin Chern number 
C sc is associated with the SHC by <j s h = C' sc = 2 in 
units of 

Now we turn on the Rashba coupling and numerically 
determine C sc . This can be done either by calculating 
each C a 'P first or carrying out the integration in Eq. 
P|) for opposite boundary condition along the i-direction 
(spin twist 9\ = — 6* = 9 X ) and common one along the 
y-direction (9 y = 9 y = 9 y ) to directly obtain C sc . We 
choose the latter method in the following. The unit cell 
of the boundary phases is divided into N mcs h = 400 to 
14400 mesh points such that the integration in Eq. is 
replaced by the sum of the solid angle Clj : C sc = • . 
Here, % = arg]l i (*jil v I / j !+ i) where i = 1 - 4 (with 
75 = ji) denote four mesh points at the j-th square of 
mesh patches in the 0-space. Numerically it is verified 
that the computed C sc is well converged and insensitive 
to mesh sizes. We find that C sc always remains quantized 
at C sc = 2 before the band gap collapses at a critical 
Rashba coupling, say, V^ = 0.33i if Vso is fixed at O.lt. 
Such a topologically nontrivial phase also survives for 
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FIG. 4: Solid angle Qj vs. spin twist 6 X and momentum k y , 
which is meshed into Nm^sh = 120 x 120 points for a pure 
system with N = 240 x 240 at Vso = O.lt and Vr = 0.3i. 
The total Berry phase ^^TJ" h ftj = 47r and thus C ac = 2. 

very small Vso- so long as Vr is not much larger than 
Vso- 

The total spin Chern number C sc can be also expressed 
as the sum of the Berry phase of each eigenstate below 
the Fermi energy. Since along the y-direction, a common 
boundary condition is imposed, which does not break 
the translational invariance of the system along this di- 
rection, one can calculate the contribution to C sc of the 
cigcnstates for a given momentum k y and formulate C 8C 
as an area integral over k y and 6 X . Interestingly, it is 
found that all eigenstates contribute more or less uni- 
formly to C sc at small Vr, which is responsible for the 
robustness of the topological invariant C sc . On the con- 
trary, at a larger Vr = 0.30£, the energy gap is about 
to collapse, and correspondingly the nonzero flj becomes 
quite localized and concentrated at four points in the 
phase space (0, ±fc*) and (it, ±fc*), as shown in Fig. 0] 
where k* = ^ (it appears that there are six peaks simply 
because of the periodicity in 6 X ). As a matter of fact, we 



find that each peak in Fig. 4 carries about half quantized 
topological charge by performing area integral around 
the sharp peak, while the rest of the area almost has 
no contribution to the Chern number C sc . With further 
increasing Vr, these half quantized topological charges 
will mix and merge with the ones coming from the upper 
energy band (with opposite signs), when the band gap 
disappears and the two energy bands touch. For a pure 
sample, a step jump of the C sc from 2 to at a critical 
Rashba coupling Vr ~ 0.33t is observed, where the en- 
ergy band exactly collapses, corresponding to a quantum 
phase transition into a metallic state. The phase bound- 
ary for the topological quantized C sc = 2 QSHE state has 
also been examined in the presence of disorder, where we 
find that the number of edge channels N e d (see Fig. 3b) 
has one to one correspondence to the spin Chern number 
C sc , and thus we identify the CNM as the topological 
origin of the nontrivial QSHE. 

To summarize, we have numerically studied several dis- 
tinctive properties of 2D honeycomb lattice model for 
the QSHE. We have identified a pair of ideal conducting 
edge states within the bulk band gap in an open system, 
which carry spin polarization and can pump spins from 
one edge to the other upon an adiabatic insertion of mag- 
netic flux into a cylindric sample. These edge states can 
be detected by using a Corbino disk of graphenc with a 
circular gatei^. Then we have established the relation of 
the QSHE and the topology of the bulk states character- 
ized by a spin Chern number C sc — ±2. The spin Chern 
number C' sc remains precisely quantized in the QSHE 
phase until the band gap collapses at strong disorder or 
Rashba coupling, where the edge states disappear simul- 
taneously. 
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